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Abstract: We present the QCD analysis of the cumulative out-of-event-plane mo- 
mentum distribution in the process pp into Zq and a hard jet (event plane defined by 
the pp and Zq momenta). Particular attention is placed on the near-to-planar events 
for which we derive the all-order resummed result to next-to-leading accuracy. We 
consider also the leading power correction originating from the fact that, even in hard 
processes, the resummed QCD coupling runs into the infrared region. We aim at the 
same level of accuracy which, in e^e~ annihilation, seems to be sufficient for making 
predictions. Contributions from a "soft underlying event" due to beam remnant inter- 
actions are discussed. Experimental data (not yet available) are needed to cast light 
on the predictive level of standard QCD analysis in hard hadron-hadron collisions. 
We plot examples of the predicted distribution at Tevatron energies. The techniques 
here developed can be extended to other hard hadron-hadron and hadron-lepton 
processes. 
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1. Introduction 

QCD radiation (jet shape distributions) in e+e^ annihilation has been intensively 
studied at the accuracy needed to make quantitative predictions [l]-[8] (double and 
single logarithmic resummations, matching with fixed order exact results and power 
suppressed corrections) . These distributions have provided various tests of QCD [9] 
and measurements of the running coupling [10]. 

For hard hadron-hadron (hh) coUisions, the analysis of the associated QCD ra- 
diation is more complex than in e^e~ annihilation. There are three main differences 
between the two processes. First, the analyzed e+e^ jet-shape distributions are 
collinear and infrared safe (CIS) quantities (all soft and collincar divergences can- 
cel). On the contrary, in hard hh-coUisions [11], jet-shape distributions arc finite only 
after factorizing collinear singular contributions from initial state radiation (giving 
rise to incoming parton distributions at the appropriate hard scale). 

A second important difference between the two processes is that, while in hard 
hh-coUisions a typical event is a multi-jet emission^ (originating from two incoming 
partons and from partons going out from the hard elementary collision), in e^e" 
annihilation the bulk of events is given by 2-jet emission (originating from the primary 
qq pair). Only recently have 3-jct event shape e+e^ obscrvables been studied: the 
thrust minor [2] and the D-parameter [3] in the near to planar 3-jet limit and 
the light-jet mass and narrow-jet broadening [12, 13] in the 2-jet limit. In this paper 
we are mainly interested in the near to planar 3-jet limit. 

Finally, in hard hh-coUisions the standard QCD description does not account for 
the entire emitted radiation. One needs to add a soft underlying event, not present 
in e+e~ annihilation, which could be considered to result from low-pt interactions 
involving the spectator partons (beam remnant interaction). The necessity of adding 
such a component to the radiation was studied in [14], in which was analysed the 
"pedestal height" in hard jet production, i.e. the mean transverse energy per unit 
rapidity accompanying a high-transverse-energy jet. The pedestal height and its jet 
transverse energy dependence measured at the CERN pp collider [15] are accounted 
for by superimposing on the standard hard QCD emission a soft underlying event 
similar to that of a minimum-bias collision. For more recent analysis on the features 
of the beam remnant interaction see [16]. 

^Actually in the particular case of processes in which only leptons are emitted in the elementary 
collision, one has just 2 jets, originating from the incoming parton radiation. 
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In this paper we consider the process 

/ii + /i2^^o + jet + ... (1.1) 

with hi, h2 the two incoming hadrons (where hih2 are pp or pp). The weak Zq boson 
is emitted with large transverse momentum Qt- Here the dots represent initial state 
jets and intra-jet hadrons. Defining the event plane as the plane formed by the hi, /i2 
and Zq momenta, we consider the distribution in the observable K^^i defined as 

i^out = E'br*h (1-2) 

h 

where p^"* is the out-of-plane momentum of the hadron h. To avoid measurements 
in the beam regions, the sum indicated by Y2'h extends over all hadrons not in the 
beam direction, i.e. emitted outside a cone around the beams. 

The process (1.1) involves three jets: the large angle jet (generated by the hard 
parton recoiling against the weak boson) and the two initial state jets (generated by 
the incoming partons). The observable i^^out is similar to T^, the 3-jet shape observ- 
able in e'^e~ annihilation, in which the event plane is defined by the thrust and the 
thrust major axes. Our analysis of Xout will then make use of the methods intro- 
duced for the study of [2]. We will obtain the following factorized perturbative 
(PT) contributions and 1/Qt power corrections: 

• incoming parton distributions at the proper hard scale. They will be obtained 
from resummations of all powers a" In^/x/ii'out (yW is the small factorization 
scale needed to subtract the coUinear singularities and, as we shall see, Kout is 
the hard scale for these distributions) ; 

• "radiation factors" characteristic of our observable. They will be obtained, for 
small Kout, from resummations of double logarithmic (DL) and single logarith- 
mic (SL) terms {a'^ln"''^^ Kout/Qt and a'^ln"' Kout/Qt respectively). Prom the 
coUinear contributions we need to exclude the pieces already accounted for by 
the reconstruction of the hard scale i^out in incoming parton distributions. As 
we shall see, the hard scales in the radiation factor, of order Qt, are determined 
by the geometry of the hard elementary process. To identify the scales to this 
accuracy we need to work at the SL level; 

• matching of the above resummed result with the fixed order exact calculations. 
This allows us to obtain a description of the distribution in the full range of 
Kout (small or of order Qt); 

• leading l/Qt-power corrections to the PT result for the radiation factor. The 
incoming parton factor has corrections which are of second order (l/Qf), see 
[5], and they will not be considered here. 
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The first three points belong to standard PT analysis. Concerning the last point, 
these non-perturbative (NP) corrections originate from the fact that, resumming 
the PT expansion, one reconstructs the running coupling as{kt) at the virtual scale 
which assumes values between kf ~ Qt and kt = 0. They should not be confused with 
soft underlying event contributions due to beam remnant interaction, which will be 
considered later. 

To deal with the running coupling in the small kf region we use the same pro- 
cedure followed in the analysis of e~^e~ jet-shape distributions. We extrapolate the 
running coupling into the large distance region using the dispersive approach [5] and 
we determine the coefficient of 1/Qt corrections in terms of a single parameter, usu- 
ally denoted by ao, which is given by the integral of the QCD coupling over the 
region of small momenta k < ni (the infrared scale //./ is conventionally chosen to 
be jii — 2GeV, but the results are independent of its specific value). Effects of 
the non-inclusiveness of Xout will be included by taking into account the Milan fac- 
tor M. introduced in [6] and analytically computed in [7]. The NP parameter cto, 
which is the same for all jet shape observables linear in the transverse momentum of 
emitted hadrons, has been measured and appears to be universal with a reasonable 
accuracy [9]. 

There are a number of differences between the analysis of -ft'out in the hh-process 
(1.1) and that of T„i in e+e" annihilation, besides that concerning initial state radi- 
ation mentioned above. First of all in the present case the event plane is defined by 
the hi,h2, Zq momenta while in e'^e" it is defined by the full structure of the emitted 
radiation. Then in e'^e~ one has a complicated interplay between the observable and 
the event plane definition (technically one has to introduce various Fourier integra- 
tion variables needed to factorize the event plane condition). These complications 
are absent in the present case. A second important difference is the presence of the 
recoil momenta of partons underlying the three jets. The recoil components enter 
the observable, the kinematics and the matrix element. In the e+e^ case all three 
primary partons generating the jets move out of the event plane, due to the recoil 
with the emitted secondary partons. In the present case, instead, the two incoming 
partons (at a low subtraction scale) are fixed along the beam direction by the kine- 
matics of the parton process. This makes relevant the presence of the rapidity cut 
excluding the beam region (see later). 

The paper is organized as follows. In section 2 we define the process under 
consideration, the observable Kout and its distribution. We specify the phase space 
region of -K^out in which we perform the QCD study. In section 3 we introduce 
the factorized structure of the distribution at parton level: the hard elementary 
partonic process, incoming parton distributions at the hard scale and radiation factor 
corresponding to our observable. In section 4 we perform the resummation at SL 
accuracy (we start from the analysis of soft contributions). We show here how to 
factorize the two contributions generating the incoming parton distributions and the 
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radiation factor. In section 5 we obtain the PT contribution to the radiation factor 
and its NP correction. In section 6 we compute the distribution, matched to the 
exact fixed-order result, and we present some numerical results. Finally, section 7 
contains a summary, discussions and conclusions. We add few technical Appendices 
A-F. In the last one (Appendix G) we discuss the contribution of beam remnant 
interaction. 

2. The process and the observable 

The incoming hadrons and Zq momenta in the process (1.1) are given by 



hi = i (1, 0, 0, 1) , /i2 = (1, 0, 0, -1) , g = (Qo, 0, Qt, Qz) , (2.1) 



with large Qt of order of the Zq mass M. The Zq is taken on-shell, and its decay 
products are excluded from our calculation. Including the hadronic decays is possible 
by following the analysis of 2-jet emission in e'^e" annihilation. The observable Xout 
is given by (1.2) where p^"* = Phx, since the event plane is the y2;-plane, and the sum 
extends over all hadrons with rapidity r)h in the range 



with 770 large. This implements a cut of angle ©0 around the two beam directions. 
To avoid a strong dependence on rjo we will consider Kgy^i not too small (see later) . 
We study the integrated distribution in Kout at fixed Qf- 



with d(Tm/dQt the differential distribution for m emitted hadrons in the process 
under consideration. We then use this to analyse the normalized Kont distribution 
for events with the Zq transverse momentum greater than some cut-off Qm'- 



with Qm a fixed upper limit. We will choose Qm at the kinematical boundary. 
3. Parton process 

At parton level, the process (1.1) is described by two incoming partons of momenta 
Pi,P2 (inside the hadrons /ii, /12), the outgoing Zq and an outgoing hard parton 
accompanied by an ensemble of secondary partons ki 



Vhl < 770 - - In tan — , 



(2.2) 




(2.3) 




P1P2 QPsh ■■■k, 



(3.1) 
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There are three configurations of the incoming partons, with pi,p2 corresponding to 
qq, qg and gq. Taking a small subtraction scale fi (smaller than any other scale in 
the problem), we assume that pi,P2 (and the spectators) are parallel to the incoming 
hadrons, 

Pi^ ^^/s{xl,0,0,Xl) , p2 = 2V^(^2,0,0, -X2) . (3.2) 
Therefore, the observable we study is 

Kout \P3x\ + l^^^l ■ (^-3) 

i 

The hard partonpa, recoiling against the weak boson, is emitted at a large angle and 
near the event plane. For small /Tout the secondary parton momenta ki are also near 
the event plane. 

The QCD calculation of the distribution (2.3) is based on the factorization of 
parton processes: 



• incoming parton distributions at the scale 



• elementary hard process; 



• evolution of the incoming parton distributions from jj, to the hard scale Kout 
obtained by resumming contributions of partons coUinear to pi and P2; 

• radiation factor, which for small K^ut is obtained by resumming contributions 
of partons soft and/or coUinear to the three hard partons; 

• soft underlying event due to beam remnant interaction. 

Concerning the last point, one expects a contribution to K^ut from the beam remnant 
interaction [14] which could be estimated by 

^remnant ^ 2r^o^f = (^^)remnant ^ ^3 ^^, 

TT 

with Af and (^}^^y^^^^^^ ^]-^g mean number per unit rapidity and the mean kt of hadrons 
produced in the beam remnant interaction (see Appendix G). From the study in 
[14, 16, 17] one estimates i^^^mnant ^^le order of few GeV. The models considered 
for the beam remnant interaction (in the central rapidity region) do not depend on 
hard scales and should therefore be the same in all hard processes. 
In the next sections we discuss only the hard QCD pieces. 
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3.1 Incoming pcirton distributions at the scale 

We denote by q[{x), q[{x) and ga{x), the distributions in the momentum fraction x of 
quark, antiquark and gluon inside the hadron ha (with a = 1, 2) at the subtraction 
scale n- The quark and antiquark carry the flavour index /. We introduce the 
incoming parton distributions at the scale fi for the three configurations 

Vlq{xi,X2,n) = q{{xi)qi{x2) + q{{xi)qi{x2) , 

p/(,x„.„;.) = [gf(.0+g-f(-i)]^2(x2), (3.5) 

Vlg{xi,X2,i^) = gi{xi) qi{x2) + 4(^2) 

3.2 Elementary hcird process 

We neglect the secondary emitted partons ki and the hard process (3.1) reduces to 

PlP2^QP3, (3.6) 

where the incoming partons, the weak boson and the outgoing parton momenta can 
be written as 

Pi = ^VSe^il, 0, 0, 1) , P2 = ^V§e-\1, 0, 0, -1) , 

Q = ( -E" cosh y+p cos ^sinhy, 0, psin^^, pcosOcoshy+Esmhy) , 

P3 = {pcoshy—pcosOsmhy, 0,—psm9, — pcos^ cosh^z+psinhT/ ) , (3-7) 

p^—j^, , Qt^psmO, 8^X1X28, y = iln— . 

2vs 2vs ^2 



Here E,p and 9 arc the energy, momentum and angle of Q with Pi in the centre of 
mass system of process (3.6); y is the boost along the z-axis bringing the process 
(3.6) from the centre of mass to the laboratory system; Xi,X2 are the momentum 
fractions of Pi, P2 entering the elementary coUision. 

The elementary cross sections dag{s,Qt)/dQt for the three partonic configura- 
tions are given in Appendix A. In the elementary process (3.6) the outgoing momen- 
tum P3 is in the event plane so i^^out — 0. Summing over flavours / and incoming 
parton configurations S — qq, qg, gq, we have the Born contribution to (2.3) as 

3.3 Distribution at parton level 

Considering the secondary emitted partons ki, the hard momentum P3 moves out of 
the event plane, acquires a soft recoil qs and the observable is 

Kout ^ ^ \kix\ + \q3x\ , Ps^Ps + qs, q3x^ - ^kix ■ (3.9) 
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For i^out ~ Qt the distribution is obtained from exact fixed order results. In the 
following we consider the region of small -K^out in which one needs all order QCD 
resummations. In this region the distribution (2.3) can be factorized as follows 

with the lower bound sq given in (A. 9). The distribution includes the incoming 
parton distribution V{n) in (3.5) and resums the (factorized) coUinear and soft pow- 
ers of ln///i^^out or lni<'out/<5t) with /i the subtraction scale for coUinear singularities. 

In (3.10) we have factored out the elementary parton cross sections da/dQt given 
in (A. 8). The coefficient Cs{as) = 1 + O (ag) is a non-logarithmic function which 
takes into account hard corrections not included in Xj. It can be computed from the 
exact fixed order results. As one expects, and as will be discussed in detail in the 
next section, the distribution Xj can be factorized as follows 

X/(s, Qt, K,,,) = J dX.dX, S (^^ - X,X2^ Vl (X1X2, i^out) • As (s, y, Qt, K,^,) , 

(3.11) 

with y — ^\nXi/X2. We have two factors: 

• the first, Vg , is the incoming parton probability evolved from /i to the hard 
scale Kout- It resums singular terms coming from secondary partons which 
are coUinear to the incoming partons pi and p2, giving rise to the anomalous 
dimensions; 

• the second. As, is the radiation factor corresponding to the observable i^^out- 
It resums powers of InXout/Qt a-nd is a CIS quantity. It is sensitive only to 
QCD radiation and therefore does not depend on the flavour (we neglect quark 
masses). There are various hard scales in As (given in terms of the invariants 
{PaPb) ~ Qt) which will be determined by the SL accuracy analysis. 

In the following we will obtain Xj by resumming the QCD radiation in the rapidity 
region (2.2). To simplify the analysis we will consider Xout sufficient large, in the 
region 

K,nt>K^t-^Me-^°. (3.12) 

As we will discuss in detail, here the PT results at SL level do not depend on rjo, 
(hadrons h' emitted inside the beam cones typically have \ph'x\ smaller than K^^^). 
The dependence on rjo and on the boost y enters only in the NP corrections of As 
(NP corrections affect the distribution at any rapidity). 

Ideally, in order for our results to be valid over the widest possible range of -fCout, 
we would like to take 770 as large as is experimentally possible. 
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4. Resummation and factorization of 

In this section we derive the factorization structure (3.11). In particular we show 
that, to our accuracy, in Vg the hard scale is actually given by Kout and the rapidity 
cut 770 is irrelevant. We also deduce the expression of the radiation factor As to 
next-to-leading order which will be discussed in the next sections. 

We consider first the resummation of logarithmic terms coming from soft sec- 
ondary partons. They include all DL terms and the SL terms originating from soft 
partons. The remaining SL contributions (collincar non-soft secondary partons) will 
be included later. They give the non-soft part of the anomalous dimension and 
contribute to fix the hard scales in As, see [2]. 

We resum the enhanced soft contributions to next-to-leading order by using the 
factorization of soft radiation. To this end we extend to the process (3.1) the methods 
previously introduced in e'^e" to analyse the distribution of CIS observables for 3-jet 
events, see [2] . The new fact in the present case is that the distribution is not a CIS 
observable and depends on the subtraction scale fi. 

The square amplitude for the emission of n soft partons in the process (3.1) can 
be factorized as follows 

\Mlniki ...kn)\'^ \mIJ' ■ Ss,niki ■ ■ ■ A;„) . (4.1) 

The first factor is the Born square amplitude which gives rise to the Born distribution 
ddg/dQt in (3.10). The second factor is the distribution in the soft partons emitted 
from the system of the three hard partons pi,p2 and P3. It depends on the colour 
charges of the emitters in the various configurations 5. 

By using (4.1), the soft contributions to the distribution Ig are resummed by 



dxidx2Vl{xi,X2,ij) XI ~i / n Ss,n{ki . . .kn) 



(4.2) 



Here the momentum fractions Xi, X2 of the hard elementary process (3.6) are given 

by 

X, = xallzi, a = 1,2, (4.3) 

where Zi is the splitting fraction associated with collinear radiation from the incoming 
partons, and is the region in which ki is collinear to Pa- As we shall discuss later, 
the precise form of the coUinear regions is not crucial. The soft factor Ss, n depends 
on the hard elementary collision variables s, y, Qt in (3.7) and the recoil momentum 
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To resum the expansion we write the constraint on s in the form 

<^-■■■^■^^i■S/Sfe)"■-■^.■■-■ .(;-.«.). 

(4.4) 

and the phase space of (4.2) in terms of Fourier and MeUin transforms 



€> ^Kout - ^ \kix\ - Ifel j S{q3x + ^ ki^) 



(4.5) 



The Melhn Na- and //-contours run parallel to the imaginary axis with Re A^^ > 1 and 
Rei/>0. Using (4.4) and (4.5), X(Xout) can be written in the form 



(4.6) 



where the distribution Ss{i^, Ni, N2, /i) is obtained by resummming the soft contri- 
butions 



n ' i=l ^ * 



Here the source U{k) takes into account the phase space constraints (e'^^'^'^) and the 
observable (e~''l'^"^l) for rapidity rjk in the region (2.2) 

U{k)=uik^) = e-'^d'^^l-*/^'^^) for|%|<77o, 
U{k) = Uo{k,) = e'^^'^ for |%| > 770 ■ 

The source 612(2;) takes into account the fact that the coUinear radiation reduces the 
longitudinal momentum components of the incoming partons pi and p2- 

ei2{z) = 2;^"-^ , for /c e C, a = 1, 2 ^^^^ 
euiz) = 1 , for /c ^ Ci,C2 ■ 

The near-to-planar region K^ut <S Qt corresponds to the region of the MeUin variable 
uQt ^ 1. Having introduced the Fourier variable ufi conjugate to k^, the condition 
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Nc 
2 




Nc 
2 




Nc 
2 



\kx\ ~ -^out ^ corresponds to /5 = (in other words, the /3-integration will 

be fastly convergent). Therefore we have here that pQt is the only large parameter 
we need to consider. 

To obtain the exponent TZs at SL accuracy we follow the same procedure de- 
scribed in detail in [2] and for the configuration 5 we have 

Ws(k) [1-U(k)eu(z)]. (4.10) 

Here Ws{k) is the distribution of soft gluon radiation off the hard three-parton an- 
tenna in the configuration S given by 

1 

Wl3 + W23 - J^W 

Wi2 + ^23 - -^Wi3 ) , (4.11) 
1 

where Wabik) the standard soft distribution for the emission of a soft gluon k from 
the a6-dipole 

Here the running coupling is taken in the physical scheme [18] and kab,t is the invari- 
ant transverse momentum of k with respect to the Pa,Pb hard partons. The unity in 
the square bracket in (4.10) takes into account the virtual corrections. The expres- 
sion in (4.10) resums all enhanced terms at ncxt-to-leading order coming from soft 
contributions. To reach the full SL accuracy, one needs to take into account also the 
non-soft part of the coUinear splitting which will be considered later. 

The exponent TZs in (4.10) differs from the e+e~ radiator in 3-jet events by the 
presence of the z^""^ factor in the sources ei2 when k is coUinear to one of the 
incoming partons pi or p2- As a consequence, TZg is not a CIS quantity and then 
depends on the subtraction scale 

4.1 Factorizing incoming parton distributions and radiation factor 

In the present formulation, the factorization structure (3.11) results by splitting the 
source as follows 

[1 - U{k)e^2{z)\ = [1 - U{k)\ + [1 - e^2{z)\U{k) , (4.13) 
so that the exponent TZ can be split into two terms 

7^^(l/, TV,, ;,) = Rs{u, (3) + Ts{v, TV,, N^, /x) . (4.14) 
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The first term, which produces the radiation factor in (3.11) is given by 

— Ws{k) [1 - U{k)] . (4.15) 

It is a CIS quantity, independent of N^, of the same type as the radiator in 3-jet 
e'^e~ processes [2, 3]. It depends also on the hard variables s,y,Qt in (3.7), on the 
rapidity cut rjo and on the recoil component q^x- 
The second term is given by 

^ r d^k 

Ts{v, /5, iVi, N2,ii) = J2 — Ws{k) [1 - z""^-'] U{k) . (4.16) 

a=l '""'>■ 

The integration is confined to k within the region Ca coUinear to Pa, so that we can 
neglect the dependence on q-^x-, in the soft limit. is coUinear singular and we 
therefore need to introduce here the subtraction scale This term gives the (soft 
part of the) anomalous dimensions of the two incoming partons, and so evolves the 
incoming parton distribution 'P(/u) to ViKout) in (3.11). 
We discuss first and then the radiator R^. 

4.2 Incoming pcirton evolution at the hard scale 

We first observe that, to our accuracy, F^ does not depend on the rapidity cut (2.2). 
This is shown by the following argument. For large vQt we have i/ ~ ~ i^Z^x so 
that the difference between the two sources 

u^{kx) - u{kx) = e^'f"'^ (1 - e-'^l^^l) , (4.17) 

vanishes unless ~ i^^out- Therefore, in the integral giving F^, we can replace the 
source U{k) with u{kx) (or with Uo{kx)) with correction of order as(-f^out)- The result 
is then independent of the rapidity cutoff rjQ, to our accuracy. 

The calculation of F^ is rather standard and is performed in Appendix B. The 
crucial point here is the identification of the hard scale. In our case the scale is 
obtained by using the fact that, within next-to-leading order accuracy, the /ca;-source 
can be replaced by an effective cutoff (see for instance Appendix C of Ref. [3]) 

[1 - u{kx)] ^ Q{\kx\-p-') , p = e^^i/^TT^ . (4.18) 

Prom this one has that the hard scale is of order p"^. Prom Appendix B one has 

Ts{v, ATi, N„p)c^-J — js{Ni, N2, as{k)) , (4.19) 



Here 75 is the soft part of the anomalous dimensions of the two incoming partons 

2 

^s{N^,N2,a,ik)) = J^^i'^^iiNa) , (4.20) 

TT ^ — ' 

a=l 
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where the function Si{N) is given in (B.2) and 

Cf^Cf'^Cf^Cl' = Cl' = Cl'' = CF, Cf" = CI' = Cf = Ca, (4.21) 

are the colour charges of the parton #a in the configuration 5. The soft part of 
the anomalous dimension (4.20) is accurate at two loop order provided we use the 
coupling in the physical scheme [18]. It is diagonal in the configuration index since 
soft radiation is universal and does not change the nature of incoming parton. 
Upon integration over Xa and the Mellin variables Na one obtains 

(4.22) 

The dependence cancels in the product V{fi) -e^^^^^ giving the distribution V{p^^) 
with p given in (4.18). For small K^nt "we have ~ -f^out and wc can replace 
the hard scale with i^Tout with corrections of order as(-^out) not enhanced by 
logarithms. 

Up to now we have considered only soft contributions giving (4.20), the leading 
part of the anomalous dimension for large Na- However for the process under consid- 
eration we have to consider contributions for — O {1) {Xa not too close to one) 
and we need to consider the full anomalous dimension which is no longer diagonal 
in the configuration index 5. In this way the incoming partons may change from a 
quark to a gluon and vice versa. The resulting exponent F becomes a matrix in 6 so 
that the product V{p) ■ e'^^/^) in (4.22) becomes a matrix product and one obtains 
the parton distribution {Xi, X2, Kout) fully evolved from the subtraction to the 
hard scale. We do not consider here power corrections since they are of second order 
(1/Q?), see [5]. 

Finally, (4.22) gives the factorized result (3.11). 



5. Radiation factor 

The distribution As{s,y,Qt, Kout) entering into the factorized expression (3.11) is 
obtained from the piece Rs{u,P) of the full radiator in (4.10). From (4.7) we obtain 

-^-T^ As{s,y,Qt,i^) , (5.1) 

where 

Ms, y, Qt, = r dqs, r ^ e-'^d^^^l-^'^'^-) e-^^^'^'^) . (5.2) 
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We recall that the CIS radiator Rs{iy, 13) depends on the elementary collision variables 
s, Qt (see (3.7)), rapidity cut rj^ and recoil component q^x- The radiator Rs contains 
a PT contribution and an NP correction 

Rs{v,(3)^R^,'^{y,(5) + 5R5{v), (5.3) 

which we discuss in the following. 

5.1 The PT radiator 

Here we obtain R^'^ to SL accuracy. We start from (4.15). The contribution from 
the hard coUinear splitting will be considered later. 

First of all, we can neglect the soft recoil ^3 since, as shown in [8] , it contributes 
beyond SL accuracy. Then we set Pa Pa (see (3.7)). Moreover, to SL accuracy, the 
PT radiator is independent of r^o and as long as one considers i^out in the region 
(3.12). This is due to the fact that for \rik\ > f]o we can replace Uo{k) with u{kx)- 
The argument is similar to the one presented for the evaluation of (see Appendix 
C for a detailed discussion on this point). 

The calculation of the PT radiator is performed in Appendix C and, to SL 
accuracy, one finds 

= t (p, Ql) + Cf t (p, QJ-) + Ca t (p, Q^) . (5.4) 

where 

( n^- dk2aj2k) r-^ 
r {P, Qa) = J^ ^ y In - , Kr.^^, (5.5) 

with as in the physical scheme [18]. The factor 2 in the argument of the running cou- 
pling takes into account the fact that the observable involves only the x-component 
of transverse momentum, while the ^/-component is integrated out. The upper limit 
K in the k integration is beyond our accuracy as long as of order of the hard scale 
\/^. The PT hard scales are given by 

Qf = Qf = Q12 ■ e-'/' , 

Qf = Qf = Qis ■ e-'/' , 

= = Q23 ■ e-'/' , 

The first factors, expressed in terms of the Qab in equation (C.4), are determined by 
the large-angle soft emission. The quark or antiquark scale is given by the invariant 
mass of the quark- antiquark dipole. The scale for the gluon is its transverse mo- 
mentum with respect to the quark-antiquark dipole. The rescaling constants e~^/^ 



.gg _ Ql3 -Po/iNc 



= :il±^ . e-^o/*^^= , (5.6) 



.gq _ Q12 Ql3 -/3o/4JV, 
V23 
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and e~^o/^^'= take into account SL corrections coming from the hard parton sphtting 
functions. These constants and the precise expression of the geometry dependent 
scales Qa are important only beyond DL accuracy. In conclusion, the PT radiator 
depends on the geometry of the event, on the underlying configuration 5, and on the 
presence of the recoil in the kinematics (/3-dependence) . It does not depend on the 
boost y since the rapidity cut r/o does not affect the PT result. 

5.2 NP corrections to the radiator 

The procedure for computing the leading NP corrections, including two loop order 
to take into account the non-inclusiveness of jet observables, is the usual one [6] and 
recalled in Appendix D. One finds 

5i?,(g3.) 1^ (c? (7?o - r/3) + Ci'^ ivo + Vs) + C? In f^) , C = 26"^ , 

V mx\/ 

(5.7) 

where A^^ is the NP parameter given in (D.6). It is expressed in terms of the integral 
of the running coupling over the infrared region 

(^oilJ'i) = / — o^s{k) . (5.8) 

This parameter q;o(/^/) is the same as enters the e+e~ jet shape variables 1 — T, C, 
M'^/Q^, B, Tm and D, see [2, 3, 4]. After merging PT and NP contributions to the 
observable in a renormalon free manner, one has that the distribution is independent 
of Hi. 

The quantity 773 is the rapidity of the outgoing hard parton #3 in the laboratory 
system (3.7) 

^3 = lln^^ = y + ln^ = y + lntan^ . (5.9) 

Due to the symmetry of the integrand of (3.10) for collisions of interest to us, we may 
take cos^ > as discussed before (A. 8). Finally, the colour charge Ca^ of parton 
#a in the configuration 5 is given in (4.21). 

The result has a simple interpretation based on the fact that the observable is 
independent of the soft gluon rapidity. First consider the contribution from the out- 
going parton #3, proportional to Cf \ Here the In \qzx\ contribution comes from the 
rapidity integration along the outgoing parton direction. Real-virtual cancellation, 
which takes place when the angle of the outgoing parton with the event plane 
exceeds the corresponding angle of the soft gluon, provides an effective rapidity cut 
and leads to the In |g3a;| term. The rescahng factor C, is due to the fact that rapidity 
is related to the angle between two vectors while the boundary here is given in terms 
of the angle between a vector and a plane. 
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Then consider the NP correction due to emission from the incoming parton 
proportional to Cf^ and . Since the pi and p2 momenta are fixed, no recoil is 
present and the soft gluon rapidity 77^ is bounded by r/o- The r/fc-integration gives 
then (770 — 773) and (770 + 773), for a = 1 and a = 2 respectively, i.e, the length of the 
rapidity interval between 773 and the boundary ±7/0 • If we remove the 770-bound (i.e. 
we keep in the observable all partons including the ones in the beam direction) the 
rapidity of the soft gluon can go up to the kinematical limit and one obtains NP 
corrections involving the logarithmic moment of the running coupling in the infrared 
region. 



6. Distribution 

We are now in the position to obtain the full distribution E((5m, -f^out) in (2.4) to the 
standard QCD accuracy. First we obtain the resummed PT expression T^^ . Then, 
using the exact result of the O {al) matrix element calculation, we compute the 
first correction of the coefficient function C{as) in (3.10) and perform the matching 
of the resummed and the exact result to this order. We then include the leading 
power correction coming from the NP part of the radiator (5.7). Finally we plot 
this distribution for pp collisions at the Tevatron (the contribution of the underlying 
event due to the beam remnant interaction can be taken into account successively 
as a "rigid shift" in the argument by the quantity i^^Jj^™"^"* in (3.4)). 

6.1 Resummed PT contribution 

The PT contribution TF^ [Q.^^ Koat) to SL accuracy is obtained from the radiation 
factor (5.1) by taking only the PT part of the radiator given by (5.4). Performing 
the Mellin transform (see Appendix E) we obtain, to SL accuracy, 

^r(S,Q„ifou.) = e-™.)-^|^. (6.1) 
Here i?^^ (-f^out) is given in (5.4) with p replaced by K~^^ and 

TT Aout (1 + ^R') 

(6.2) 

with K the same hard scale as in (5.5) so that 

Ar{s,Qt,K) = l. (6.3) 

To first order in ctg we have 

ATis, Qu Km) = 1 - ^I^TLi + 2 In ^ + O {al) , ^ In 

(6.4) 
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The PT contribution to the (normahzed and integrated) distribution in (2.4) is given, 
to SL accuracy, by 



1 i-Qm r [ 

xC5(as) (s,Qt,is:out)}, 



out ) 

(6.5) 



with 



<T{Qm) = dQt j dXidXsJ^ |^P/(Xi,X2,X)| . (6.6) 

The exact value of the hard scale of V in (6.6) is not important, as long as of order 
-\/^; a variation can be absorbed into a modification of the coefficient Cs- In order 
to simplify the calculation of Cs for the exact fixed order expression of E in (6.6) we 
have fixed this scale at the same value K \n (5.5) at which the radiation factor (6.3) 
is one. 

6.2 Matching resummed with fixed order results 

In the matching procedure one starts by determining the coefficients of 

C,(as) = 1 + ^ci + . • • (6.7) 
from the fixed order exact results. Since only the first loop term 

j^exax^t ^ 1 + ^ ^* + . . . (6.8) 

is known we can determine only c\. The first term E^'''^'^*, which has the DL and SL 
structure 

T.T^\Qm. i^out) = -2Ct + G'n(g^) L + Ci , L = In ^ , (6.9) 

-f^out 

is obtained^ by using the numerical program DYRAD [19]. First of all we check that 
the DL and SL terms coincide with the result of our calculation, see (F.5). Then we 
compute ci which is given by 

CliQm, i^out) = -TliQm) + {^T'^'iQm, i^out) + 2Ct - Gn{Qm) L} , (6.10) 

with ri coming from the mismatch between Lk and L, see (F.2) and (F.5). The 
logarithmic terms are completely subtracted and Ci is finite for K^ut 0. This is 
shown in Fig. 1 for some values of Qm- For simplicity here we set the hard scale K 
at M so that Lk — ^ L and ri{Qm) 0. 

^Since K in the radiator (5.4) may depend on the integration variables, the logarithmic variable 
for the full distribution S is chosen to be L and not Lx- We need then to take into account the 
mismatch between L and L^- 
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Figure 1: The first order coefficient function ci as a function of i^^out for the two values 
Qm = 50, 100 GeV and K = M. 



Given ci in (6.10) we now reconstruct the matched result. For instance, using 
the so called "log-R matching prescription" , to first order we have 



mat 



(6.11) 



It is straightforward to check that the distribution TF^ obtained by using this coef- 
ficient function reproduces the exact result in (6.8) and accounts for all terms of the 
form a"L^"~^ in the resummed expression (6.5). In order to obtain also the a^L^"~^ 
terms one should perform a second order matching, which requires the fixed order 
exact result to order a^, not yet available. 

Actually this expression for the distribution S is not yet normalized to one at the 
maximum value of -ft'out- Indeed we used a normalization point at ii'out = K, 
see (6.3). The standard way to achieve the correct normalization is to substitute 



K 



K 



K 



K 



+ 1 



L = In^ 



M 



(6.12) 



in all previous expressions (except in T,f^'^^{Qm,Ko„t) in (6.10) which is already 
correctly normalized). The variable -ft'out goes to the correct kinematical boundary 



out 



K^i) for K, 



out 



K and tends to K^ut for small values. 



6.3 Including the NP correction 

We now discuss the full distribution including the leading NP corrections coming 
from the radiation factor. Recall that power corrections from V are subleading. The 
analysis is similar to the one in [2]. The NP radiator 6Rs is proportional to the 
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Mellin variable see (5.7), and so it produces a shift of /Tout 



Xout ^ K^t = Kont - SK,,t SK,,t = A5(Xo,0 . (6.13) 

The final expression, including NP corrections, is then obtained from the PT result 
of previous subsections in which we replace K^ut with K'^^^^. 

Here A^^ is expressed in terms of the NP parameter ao{ni), see (D.6). To 
evaluate As we observe that (5.7) contains a In |g3a.| term. The qs^ distribution is 
given by the radiation factor (5.2) which leads to \qsx\ ~ -f^out, so that a In Igs^l term 
produces a In Kout contribution. We find (see Appendix E.2) 



As{K,^,) = CPirjo - m) + Ci'\m + %) + In ■ 

(6.14) 



Kout 



ing^p = In eg, + ^E + i^{i + R') + (i + f) - (^) ' 

with K given in (6.2). Notice that, expanding in powers of i?', the pa contribution is 

ln^ = ln^ + C»(i?') . (6.15) 

The factor 2 here is simply due to the fact that ps acquires a recoil which is equal to 
/s:out/2 for small B! . 

The effect of the substitution in (6.13) is a deformation of the PT distribution. 
First of all the quantity (6.14) depends logarithmically on Xout (both explicitly and 
through the SL function i?') . This imphes that the PT curve is shifted by an amount 
which decreases with increasing Xout- Moreover, A^iKout) depends also on the ra- 
pidity y distributed according to the incoming parton distributions Vs[Kou\)- 

Here the situation is different from the case of broadening [8] or [2] in e+e" 
annihilation in which one obtains very singular contributions to the shift of order 
1 / The difference is due to the different kinematical situations. In the two e+e" 
cases one has to consider contributions in which some hard partons are forced to stay 
in the event plane. Then the PT distribution is given by a Sudakov form factor and 
then the 1 / contribution comes from the integration of the logarithmic term in 
the hard parton recoil. For the present observable instead, pa is not kinematically 
forced into the event plane and then its PT radiation factor reproduces a logarithmic 
contribution. 

6.4 Numerical analysis 

We report here some numerical results. We consider pp collisions at Tevatron energy 
(-\/s = l.STeV) for some typical values of Q^. Data on the i^out distribution are 
not yet available. The results depend on the two parameters Q;^jg(M^) and Q;o(/^/) 
(with — 2 GeV) which values we fix in the range determined by the 2-jet shape 
analysis [20]. 
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Figure 2: The differential distribution in (6.16) as function of -fCout- We plot both the PT 
resummed and matched result and the full distribution including NP correction for Qm = 
50GeV. Here we have taken a^jg(Mz) = 0.120 and the NP parameter ao(2GeV) = 0.52. 
The contribution due to the beam remnant interaction is not included, it is given by a rigid 



for two values of Qm- The PT curve represents (6.5) with the coefficient function 
given by (6.11). We have performed the substitution (6.12) in order to take into 
account the correct normahzation at the kinematical boundary. The rapidity cut is 
set at rjQ = 2. For simplicity we have taken the hard scale K at M. We used the 
incoming parton distribution of [21]. 

The PT+NP curve is given by the above PT expression in which we make the 
substitution K^nt into K'^^^ according to (6.13). The leading NP correction is deter- 
mined by the single parameter ao{fii). 

As discussed before, the effect of the NP substitution (6.13) is a deformation of 
the PT distribution. In particular the PT peak is shifted by about 10 GeV. This 
effect, due to the QCD running coupling in the infrared region, has to be contrasted 
with the contribution from the beam remnant interaction which corresponds to a 
"rigid shift" of the hard QCD result by an amount Kl'^'^''^'^^ of order of few GeV, see 
[14, 16, 17], proportional to 770 but independent of the hard scales (see (3.4)). 

7. Discussion and conclusion 

The aim of the present study is the understanding of the structure of radiation in hard 



shift in Kant by i^^X""""* of the order of few GeV, see (3.4). 



In figs. 2, 3 we plot the differential distributions 

d^iQm, Kant) 



(6.16) 
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hh-collisions within the standard QCD treatment. We have introduced the jet-shape 
observable Kont which is the extension to hh-colhsions of in e^e~ annihilation. 
The hard QCD analysis includes next-to-leading order PT resummation, matching 
with exact first order result, and leading NP power corrections (arising from the 
fact that the running coupling argument runs into the infrared region). To avoid 
measurements inside the beam direction we have included a rapidity cut r^o (see 
(2.2)). For values of Kout much smaller than K'^^^ (of the order of 10 GeV for 
TjQ = 2) there are powers of In i^out /-K^^ut- calculation we consider the region 

Kout > -K^out so that we do not need to resum them. As a result the PT contribution 
does not depend on rjo. The dependence on rjo enters only in the NP correction which 
corresponds to the substitution (6.13) in the argument of the PT distribution. 

We discuss now some of the features of our hard QCD result. 

In the present calculation we have various hard scales. It is then important to 
identify the specific hard scales in the various factorized pieces of the result (6.5). We 
find that Kout is the hard scale for the incoming partons distribution Vs while it is 
the lower bound for frequencies contributing to the radiation factor As- This result 
is based on the different role of real/virtual cancellations in the coUinear singular 
quantity Vs and in the CIS radiation factor As- Technically, in the present treatment, 
the factorization {Vs ■ As) results from (4.13) while real/virtual cancellation from 
(4.18). The hard scales for As are given the scales of the elementary 

hard process. They are identified at SL level and depend on the geometry of the 
event (s- and Q^dependence) and on the configuration 6 = qq,qg,gq for the two 
incoming partons. Therefore, the shape in Ko^t of the distribution (at fixed s and 
Qm) depends on the weight of the various configurations. By changing s and Qm one 
may be able to study the three configurations separately. 
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Leading power corrections come from the radiation factor (power corrections in 
the incoming parton distributions are subleading [5] and were not considered). They 
enter as a shift in the PT distribution given by (6.13). This is a feature common to all 
observables linear in the transverse momentum. The strength of the power correction 
is given by the NP parameter A'^^ (expressed in terms of aodij), see (D.6)), the 
same as introduced for the e+e~ jet-shape observables. The structure of the shift 
is characteristic of the fact that the observable Kont is uniform in rapidity, see the 
discussions after (5.9) and (6.15). The NP shift is larger then the corresponding NP 
shift for 2-jet observables since it takes contributions from three hard partons, one 
of which is a gluon. This was also the case for the near to planar 3-jet observables in 
e~^e~ [2, 3]. We then expect that higher order NP effects may come into play near the 
peak of the distribution. This calls for a deeper analysis that would address higher 
power corrections, for example, along the lines of Korchemsky-Sterman approach 
which was recently developed for some 2-jet observables in [22]. The comparison 
with experimental data (not yet available) would shed light on this important point. 

Numerical programs for exact results on the matrix elements for the process (1.1) 
are available [19] to order a'^. Thus we have been able to compute only the first term 
of the coefficient function Cs{as)- 

In e+e^ annihilation the hard QCD analysis described above has been shown 
to be sufficient to make quantitative predictions and study the universality of NP 
effects. In hh-coUisions however we also need to take into account contributions 
coming from the soft underlying event due to beam remnant interaction. With the 
hypothesis that such contributions are factorized and independent of the hard scales 
(see [14, 16, 17] and Appendix G), they give rise to a "rigid shift" in Kout by the 
amount K^^^^^* (see (3.4)). This quantity, of the order of few GeV, is proportional 
to 7]q and the minimum bias parameters M and (/c^)'^'^™"'^"*. The last two parameters 
should be the same in all hard hh-processes and could be determined and checked in 
the study of various observables. 

This completes the analysis of our jet-shape observable in hard hh-coUisions. 
The results depend on the two parameters as{Mz) and ao{ni = 2 GeV). The impor- 
tant question is then whether the present QCD standard treatment is sufficient to 
reproduce the data in hard hh-coUisions. In particular, our results should provide 
solid ground to study whether there is any need for non-hard QCD contributions. 
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A. Elementary cross sections 

We present the three elementary cross sections relating to the configurations qq, qg 
and gq in terms of the kinematical variables s and 9 introduced in equation (3.7). 
The parton-level cross section ctJ for quark fiavour / and configuration 5 is 

H = Ysj ^Lips[PiP2 ^ QPa] • |m/|2 , (A.l) 

where the Lorentz-invariant integration measure, after eliminating a trivial azimuthal 
angle, becomes simply 

/ dLips[PiP2 ^ QPa] = dcos9. (A.2) 

For the matrix elements M/ only the simple 2^2 tree-level diagrams are 
required. After averaging over colours and spins we obtain 



^ (PiP.)(P.P3) ^^-'^ 

f ,2 _ 327r^asa;T^ {P^Qf + (PsQ)^ 

' (PlP2)(PiP3) ' 

where a/ is the electroweak coupling of the quark to the Zq 

Using (3.7) to express these matrix elements in terms of s and 9 quickly yields 
the differential cross sections 



(A.5) 



dcos9 Nc ps^/^sin^6' 

d&(g{s, 9) ^ 2na,afTR {E - p cos 9f + Ap'' 

dcos9 ~ Nc ' s^{l-cos9) 

ddlq{s,9) _ 27rasafTR {E + p cos 9f + Ap^ 

dcos9 ~ Nc ¥(T+~cos9) 

We however are interested in d&g/dQt, which is given by 

— °-— =/ (icos6'()((5t -psm6') — — . (A.6) 

dQt J-i d cos 9 
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On integrating over 9 there are contributions from cos e = ±Vl - <9tVp^- Thus 

da^{s,Qt) ta.n6 da^ 



dQt 



p d cos 9 



tan 6' da^ 



cos 0>Q 



p d cos 9 



(A.7) 



cos 6<0 



Now we make use of the fact that, for proton-proton and proton-antiproton coUisions, 
these two terms contribute equally to (3.8). (This would not of course be true for 
collisions of two arbitrary hadrons.) Therefore we may take cos^ > and write 



sin 9 = Qt/p , cos 9^ \Jl-Ql/p' ■ (A.8) 



da^{s,Qt) _ ^tan 9 da^{s, 9) 
dQt p dcos9 

The relation sin 6* = Qt/p gives a lower bound on the total centre of mass energy 



s: 



(A.9) 



Finally, a consideration of the kinematics shows that, for fixed partonic energy y/§ 
the observable is bounded by 



Kl,<s + M'-2^s{Ql + M^), 



(A.IO) 



and since Qt can be as low as Qm and s can (with small probability) be as large as 
s, we have the absolute upper bound 



Kl = s + M'-2^s{Ql + M'). 



(A.ll) 



B. Incoming parton evolution at the hard scale 

We denote by Sab the contribution to Fg from a given a6-dipole distribution Wab, see 
(4.11), and consider first the 12-dipole contribution 

«i2-y: / — ^l2(^)(l-.^"-^)^^(M=y:^lW / ^^e-^^d^-w, 



fB.ll 



with 



1- 

S,{N) ^ dz — 
Jo 



N-l ^ 
n=l 



(B.2) 



proportional to the soft piece of the anomalous dimension. Here both coUinear regions 
Ci,C2 contribute. 
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The contribution from the other two radiators Si3,S23 is similar. They take a 
contribution from a single coUinear region (the a3-dipole from the region Ca coUinear 
to Pa), SO that 

sas = S,{Na) / 0^ e-^(l'=^l-^^'=^) , a = 1,2. (B.3) 

Since the Sab are singular for 0, one needs to introduce a cutoff /i on the 

kt integral. We also use the fact that, to our required accuracy, the /Cj;-source can 
be replaced by an effective cutoff (see (4.18)). We are therefore required to integrate 
over a region in k^, fcy-spacc given by 0{p~^ — \k.j;\)Q{kt — fj.). where we choose /j, to 
be less than p~^. So, to SL accuracy, we write (for instance for the a3-dipole) 

Sa3^Si{Na) , (B.4) 

J/j, l^t 

where the remainder is beyond our required SL accuracy. 

As anticipated, since here one has < <^ Qt, the precise definition of the 
coUinear regions Ci,C2 is not important. 



C. The PT radiator 



The PT radiator is given, to SL accuracy, in terms of afo-dipole radiators 

ra,{v, P)= [ —wabik) [1 - e-d^^l-/^'^^)] , wabik) = , (C.l) 

where kab,t is the invariant transverse momentum of k with respect to the Pa, Pb hard 
partons in (3.7). For the configuration 5 = qq, for instance, we have 

R7i^, P) = Y P) + r2s{iy, /3) - ^ tM^, P)^ . (C.2) 

To evaluate the a6-dipole radiator Vabii^, we work in the centre of mass system of 
the a6-dipole. We neglect at this stage the rapidity cut (2.2): we will show at the 
end of this appendix that the difference is beyond our accuracy. Denoting by P*, P^ 
and k* the momenta in this system, we introduce the Sudakov decomposition 

P: = ^(1,0,0,1), p; = ^(1,0, 0,-1), r = aP: + /3P; + «:, (C.3) 
where Ql^ = 2{PaPb) so that (3.7) gives 

<9i2 = « , <9i3 = (1 + cos e) , Ql^ = p^^s (1 - cos ^) , (C.4) 
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and 9 is given in equation (A. 8). Here the two-dimensional vector k is the transverse 
momentum orthogonal to the a6-dipole momenta («;^ = A;^^ J. We have then 



Wab{k) = 5- , = , a> —2- . (C.5) 

Since, neglecting the recoil qs, the outgoing momentum P3 is in the y2;-plane, the 
Lorentz transformation is in the yz-plane and our observable kx remains unchanged. 
The a6-radiator has then the form 

Jo TTK^ TT J^/Q^^ a 

where the factor 2 comes because we have integrated only over the "right hemisphere" 
a > n/Qah- Integrating over a and Ky we have 

ra,{v, (3) = 2 ^ [1 - e-<\'^^\-'^^^^] ■ ^iM^l in ^ . (C.7) 

J -Qab ^ ^kx\ 

To show this we introduced Ky — t ■ \kx\ and used 

— ( \kx\Vl + ^] In ^ aJ2%\)l^S^ . (C.8) 



f 

J — c 



We extended the t-integration to infinity since it is convergent, then we integrated 
over t by expanding to second order. Corrections are beyond SL accuracy. Finally, 
using (4.18), we obtain 

Assembling the various dipole contributions and including hard coUinear splittings 
then yields, to SL accuracy, (5.4) and (5.5). 

We now show that the rapidity cut (2.2) is negligible for the PT radiator. The 
difference between the radiator with the cut imposed and that without is given by 
the dipole contributions 



Ar„,(i., /3) = - / — wab{k) e-^^^ (l - e-^l^^l) . (C.IO) 

J\m\>m 



In order to implement the rapidity cut, we express r]k the soft gluon rapidity in the 
frame (2.1) in the invariant form: 

where /ii,/i2 are the incoming hadron momenta in (2.1) and Pi,P2,y are the hard 
incoming parton momenta and rapidity in (3.7). 
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For the 12 dipole we have rjk — y + ln{aQi2/K), and thus we obtain from the 
"right" hemisphere 

Aruii^, (3) = - f ^ e-^^. (l _ e^l^'^l) , g;, = g,,e^-''o . 

(C.12) 

Here the scale for the correction is Q'-^2 ~ -^out ^ ^ ('^~^) ^"<i contribution 
is of order Q;s(ii'out) without a logarithmic enhancement. The same is found for the 
"left" hemisphere contribution. 

For the 13 dipole we obtain a similar result for the cut around the Pi direction, 
(using (D.13) to evaluate r^fc), while the cut around the P2 direction gives a tiny 
correction proportional to the size of the hole: the converse holds for the 23 dipole. 



D. NP corrections to the radiator 

We consider the NP correction 5rab to the a6-dipole radiator. In this case, as we shall 
see, we need to retain both the recoil q^^ and the rapidity cut 770 (recall that in the 
PT component they both gave subleading effects and were neglected) . We write the 
integral in the a6-dipole centre of mass variables a, (3 and k introduced in (C.3) and, 
to obtain the NP correction dvab-i we perform the following standard operations: 

• the running coupling, reconstructed by two loop emission, is represented by 
the dispersive form [5]. Then, the a6-dipole radiation Wab{k) is written in the 
afe-centre of mass system (see (C.5)) in the form 

• to take into account the emission of soft partons at two loop order [6], wc need 
to extend the source u{kx) to include the mass m of the soft system. We assume 
kx = Kcoscp — >■ V + rn^ cos 0, with (j) the azimuthal angle of k. Similarly we 
introduce the mass in the kinematical relations such as a(3 = {k^ + 'nn?)/Q\^ 
for the a6-dipole variables; 

• we take the NP part 5acf[{m) of the effective coupling. Since it has support 
only for small m, we take the leading part of the integrand for small and m. 
In particular we linearize the source U{k) 

[1 - U{k)] v^Jk^ + m2 I cos 01 e(?7o - |^fe|) ■ (D.2) 

Recall that 77^ is the rapidity of k in the laboratory system (2.1). Here we have 
neglected terms proportional to /3 since they vanish, by symmetry, upon the /3 
integration; 
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• the recoil component q^x of the outgoing parton does provide an effective cut in 
the soft gluon rapidity along the outgoing parton [2, 8]. This is due to a real- 
virtual cancellation which takes place when the angle of the outgoing parton 
Ps with the event plane exceeds the corresponding angle of the soft gluon. The 
detailed analysis of real and virtual pieces entails that the contribution from 
the observable \/l^~+rn? \ cos (f)\ in the linear expansion of the source (see (D.2)) 
has to be replaced by 

V^c^ + m'^ I cos 01 \/ k'^ + cos (f) + oi^-ix ~ (^\%x\ (D-3) 
with a the Sudakov variable in the a3-dipole centre of mass; 

• to take fully into account effects of non-inclusiveness of jet observables at two- 
loop order, we multiply the radiator by the Milan factor [6, 7] 

3 (1287r + 1287rln2-35^2)^ 5^2^ 

^ = 64^ lie. -2.; ' (°-^' 

using Uf — 2)] 

• the NP correction is finally expressed in terms of the parameter 

After merging PT and NP contributions to the observable in a renormalon free 
manner, one has that the distribution is independent of /// and one obtains 

A^P= CK„,, M^^iii |ao(/x/) - a, - /3o^ (in £ + ^ + I , (D.6) 
where 

a.^a^(Q), X^C^(^-^)-^n;, /^o = ^ - ^ ■ (D.7) 

The K factor accounts for the mismatch between the Ms and the physical 
scheme [18] and ao{fii) is given in (5.8). Here Q is the renormalisation scale 
used in the next-to-leading order PT calculation. 

The numerical coefficient c^^^^ depends on our observable i^out- For instance, 
the shift for the r = 1 — T distribution is 

da . , da^^ , ^ ^ a ^ Cr A^^ ^ ,^ 

— {r)^——{r-Ar), A,^Cf- , = 2 , p.8 

dr dr ck^^^ 

where Cp enters due to the fact the two-jet system is made of a quark-antiquark 
pair. 

We recall that these prescriptions correspond to taking into account NP corrections 
at two-loop order in the reconstruction of the (dispersive) running coupling and in the 
non-inclusive nature of the observable. We implement the rapidity cut by expressing 
r]k the soft gluon rapidity in the invariant form (C.ll). 
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D.l Dipole 12 

This procedure gives, for the 12-dipole contribution, 

Sri2 = [ dm^5a^s{rn)-A^ [ — ^ / ^Vk? + m?\ cos0| • 2/^0 

TT J dUl J K -\- TTl J -^TT 



(D.9) 



i/A^P 2r)o 



Wc used the 12-dipole centre of mass variables a, (3 and k introduced in (C.3). This 
result is found by using (C.ll) and 



(D.IO) 



P2k _ a _ (aQu)'^ 
Pik j3 + m? ' 

so that the a integration yields 

/ —e ho - In . ^ ^ 2770 , «;,m^0. (D.ll) 

The observable is uniform in rapidity and its integration gives 2770- Corrections 
coming from the presence of the recoil ^3^; can be neglected in this case. 

D.2 Dipoles 13 and 23 

We consider now the NP corrections 5ra3 to the a3-dipole radiators and again we use 
the a3-dipole centre of mass variables a, (5 and k introduced in (C.3). The situation 

is different from the previous 12-case in two respects. 

First of all we have that the 7]^ cut in the soft gluon rapidity 77^ does not affect 
the region along the hard outgoing parton ^3. There one has to take into account 
that, as in the case of broadening [8] or thrust minor [2], the recoil component q^x of 
the outgoing parton provides an effective cut in the soft gluon rapidity. If the three 
momenta in this system are given by the Sudakov decomposition 

^3* = %^(i,o,o,i), = ^(1,0, 0,-1), k*^aP; + (5P: + K, (D.12) 

we are required to use the expression given in (D.3) as our linearized source. 

The second complication for the a3 case is that the expression of |77fc| in terms of 
the variables in (D.12) is rather complex, due to the fact that the Lorentz transfor- 
mation to go from (3.7) to (D.12) involves both a yz-rotation and a boost along the 
y-axis. For our analysis, we need the expression for 77^ only for the soft gluon close 
to the incoming parton direction. 

We consider first the case of 13-dipole. For k nearly parallel to Pi we have P ^ a 
and this gives 

P2k _ fQuV K^ + w? 



Pik \Q13J {aQn) 
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The NP correction to the 13-dipole radiator is then given by 



/l3 



TT 

7r 



dm^daes{m)- 



dm? J K? + m? 



'13 



J ~ { v^K?+rn? cos (j)+aq3^ 6 (^0 + In 



(D.14) 



Again, in the region of k* emitted in the hemisphere {a < \/ k? + m^/Qia) the 
rapidity cut rjo gives the lower hmit of the a-integral. In the other region of k* 
emitted in the P3* hemisphere {a > y/ + m^/Qia) it is the recoil component ^33; 
which provides the upper limit of a. We have 



am < a < au , 
giving 

' .L. 27T 



(y.'m, — 



aM = 



13 



(D.15) 



da 
a 



^ Vk?+ vn? cos (t)-\-aq3x —oi\(lzx\^ 

\/k'^ + w? (rjo - y + In + In ^^^'f ) , C = 2e' 
V Q12 Wsxl J 



(D.16) 



Here C comes from the integration region of large rapidity of k near p^. In conclusion 
we have 



5ri3 = i/A ?7o - 1/ + In 



Q 



13 



Q 



In 



12 



Q13C 

\Q3x\ 



(D.17) 



A similar result is obtained from the last 23-dipole radiator: 



Q 



Q 



12 



Q23C 

\Q3x\ 



(D.18) 



A compilation of these three contributions then gives the result (5.7). 



E. Distribution 
E.l Evaluation of Af^ 

Here we compute A^"^ obtained from (5.1) by taking only the PT part of the radiator 
given in (5.4). Since R^^ does not depend on the recoil, in (5.2) we can freely 
integrate over q^^ to get 



7r(l + /32) 



r diy e'^'^out r 
J 2mv J_ 



Ar{s,Qt,Kut) = r^^^^ ^ , , u = e^-i^. (E.l) 
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We now perform the Mellin transform to SL accuracy. We make use of the 
operator identity 

G(^) = G(i^-i) (E.2) 

for any logarithmically varying function G. (To prove this, multiply both sides by the 
F-function operator and use the definition T{z) = dx x^~^e~^.) Thus we obtain 
the quantity A^"^ in the form 

(E.3) 



\ dmK, 



out ^ -'^'^out 

out , 



We make a logarithmic expansion of the radiator, neglecting contributions from the 
second logarithmic derivative, which are beyond SL accuracy. So we obtain 

RV (KutVi + P^) = RV {Kut) + {lE + In VT+p^) B! , 

To SL accuracy, R' is given by (6.2). 
Therefore we may write 

' r(l + fl') i_„,r(l + /J2)i+^B" ^ ' 

which leads immediately to the result (6.1). 

E.2 Including the NP correction 

The analysis is similar to the one in [2]. We report only the essential steps. Consider 

^^^^^ ^ l^d/3dfe ^_.(|,3,|_,g,3,)^-/lPTfpyi^) 

X 1 - A^M C['^ ivo - %) + Ci'^ iVo + %) + C!,'^ In ^"^^ 



k \ mx\ 

where we have expanded exp(— ^it!^) to first order in order to obtain the leading 
correction. Performing the q^^ integration we get 

^(z/)=iP»-z/A^P/5(z/), (E.7) 

where 

r°°dBp~^^^i''^^^) / 

(E.8) 
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with xiP) = In ^1 + /?2 + /?tan-i (3. This gives 

(E.9) 

Prom equation (E.2) and (E.9) we obtain 



/OO ' ^ ^ 

■CO ni 



\NP /-oo j/o tl— 



X e-<^(^outAA^^) _ + ^-2(770 + 773) + c,[HCK^M + xm) ■ 

(E.IO) 

Neglecting contributions from the second logarithmic derivative of R^^, which are 
beyond SL accuracy, we obtain 

(E.ll) 

Performing the /3 integral then gives 



_ \NP 

(E.12) 



in other words the distribution is shifted by — A^^A^, where A^ is given in (6.14). 
F. Matching 

Expanding the numerator of the integrand of (6.5) we have (flavour and configuration 
indices are understood) 



C(a,) • V{X^, X2, Xout) • A{s, Qt, Xout) = V{X,, X^, K) 

I. ZTT J -flout 

where we have expanded around -f^out = M. Here Ci is the first order term of the 
coefficient function C{as) in (6.7) and we have 

(F.2) 



' M ^ \ M K J 
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The quantity 'y{Xi,X2, K) is defined by 

-P{X„ X„ X) . 7 (F.3) 

Therefore the first order contribution Ei in (6.9) is: 

j:i^-2CtL^ + GuL + Ci, Ci = ci + ri, (F.4) 

with 

GuiQm) = (Gn) , n = (fi) , (F.5) 
where we introduced the averages 

{F) = a-\Qm) dQtJ dX^dX2Y^^-^V{{X,,X2,K)■Fsj . (F.6) 
The normahzation a{Qm) is given by (6.6). Notice that ri = for X = M. 
G. Particle production from beam remnant 

Here we discuss a factorized model for the beam remnant interaction which gives the 
contribution X^^^"^"* to the shift as in (3.4). We assume an independent emission 
model, roughly similar to the one discussed in [14, 17], in which the two outgoing 

— * 

hadron remnants produce soft particles with distribution in kt and rapidity given by 

J„ remnant Kflfl 9 

'fZ _ ll^„-6fct /T, \remnant _ ^ /p . x 

where H is the number per unit rapidity of soft hadrons emitted by the beam rem- 
nants and (j^^Y^"^'^'^'^^ ^Jig mean value. Then we obtain an additional term in the 
radiator due to the beam remnant 

^remnant I _ ^-vkt\cos4>\^^ I 

2b f 2b' , [b^ ^^-^^ 



7r(62-i/2) y^l)2_j^2 y b+U 

For i/ ~ 6, i.e. /Tout ~ ^/^^^''emnant distortious of the distribution. For u <^ b, 

i.e. Kout ^ {ktY^^'^^'^^ , we may expand 

^remnant ^ ^ ^V0^[ ^-1 ^ _ _ _ _ ^ ^remnant ^ _ _ ^Q g^ 

TT 

with X^'^^'^^^ given by (3.4). We conclude that, in the region we consider, the 
factorized beam remnant interaction can be taken into account simply as "rigid 
shift" 

i^out ^ i^out - K—' , (G.4) 
of the hard QCD distribution. 
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